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ABSTRACT 



Context. Recent asteroseismic observations have led to the determination of rotational frequency splittings for t = 1 mixed modes in 
red giants. 

Aims. We investigate how these observed splittings can constrain the modelling of the physical processes transporting angular mo- 
mentum in stellar interiors. 

Methods. We first compare models including a comprehensive treatment of shellular rotation only, with the rotational splittings ob- 
served for the red giant KIC 8366239. We then study how these asteroseismic constraints can give us information about the efficiency 
of an additional mechanism for the internal transport of angular momentum. This is done by computing rotating models of KIC 
8366239 that include a constant viscosity corresponding to this physical process, in addition to the treatment of shellular rotation. 
Results. We find that models of red giant stars including shellular rotation only predict steep rotation profiles, which are incompatible 
with the measurements of rotational splittings in the red giant KIC 8366239. Meridional circulation and shear mixing alone are found 
to produce an insufficient internal coupling so that an additional mechanism for the internal transport of angular momentum is needed 
during the post-main sequence evolution. We show that the viscosity v add corresponding to this mechanism is strongly constrained 
to be v ac id = 3 x 10 4 cm 2 s -1 thanks to the observed ratio of the splittings for modes in the wings to those at the centre of the dipole 
forests. Such a value of viscosity may suggest that the same unknown physical process is at work during the main sequence and the 
post-main sequence evolution. 

Key words, stars: rotation - stars: oscillation 



> 

m 

(N 
O 

o 



X 



1. Introduction 

Rotation is an important process that can have a significant im - 
pact on the evolution and properties of stars (e.g. Maeder 2009). 
The effects of rotation and in particular the transport of chemi- 
cals and angular momentum by meridional circulation and shear 
instability hav e been included in stel lar evolution codes (see 
the review by Maeder & Mevnet 2012). This has generally been 
done in the context of shellular rotation, which assumes that a 
strongly anisotropic turbulence leads to an essentially constant 
angular velocity on the isobars dZahnl [l992). To progress in our 
understanding of the dynamical processes at work in stellar inte- 
riors, one needs to confront these theoretical models with obser- 
vational constraints. 

Surface abundances are used as diagnostics of the efficiency 
of rotational mixing in stellar interiors. In the case of low-mass 
stars, measurements of light element abun dances are particu- 
larly valuable for this pu rpose (see e.g. iTalon & Charbonnefl 
insonneau iH l2010l) . Rotation is indeed found to have 
an impact on the lithium depletion o f a solar-type star from 
the very beginning of its evolution (Egg enberger et all 120121) 
to more advanced phases as a red giant dPalacios et al.l 12006). 
Interestingly, the predicted lithium abundance at the surface of 
main-sequence stars is very sensitive to the modelling of rota- 
tion. While lithium abundances observed for stars on the blue 
side of the lithium dip can be reproduced by rotating models in- 
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eluding shellular rotation dPalacios et al.l200"3l) . the same models 
fail at repr oducing lithium abundances for s tars on the cool side 
of the dip (ITalon & C haibon nel[l99ll2003|) . This indicates that 
there is an additional physical process operating in the internal 
layers of low-mass stars, which are efficiently spun down via 
magnetic braking during their evolution on the main sequen ce 
(ITalon & Charbonnell 19981 120031: ICharbonnel & Talonll2005l) . 

Additional constraints of the modelling of angular momen- 
tum transport in stellar interiors come from measurements of the 
surface rotation of stars. Observations of rotational periods for 
young solar-type stars in open clusters suggest that slow rota- 
tors develop a high degree of differential rotation between the 
radiative core and the convective env elope, while solid - body ro- 
tation is favoured for fast rotators (Irwin et al ■ l|2007t iBouvieri 
2008; iDenissenkov et al.ll20ldl lDenissenkovll2010l) . Models in- 
cluding a detailed treatment of shellular rotation show however 
that meridional circulation and shear mixing alone do not cou- 
ple the core and envelope, even in the presence of rapid rotation 
(Egg enberger et al1l2010al) . This also indicates the need to intro- 
duce an additional mechanism for the transport of angular mo- 
mentum during the beginning of the main-sequence evolution of 
solar- type stars. 

In addition to observations of stellar surface properties, 
those of the solar five-minute oscillations have provided 
a wealth of information on the internal structure of the 
Sun and have led to the determination of its rotation pro- 
file dBrown et alJI 19891: iKosovichev et alJ[l997l: ICouvidat et all 
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I2003t iGarcfa et alJ l2007h . This major constraint on the mod- 
elling of angular momentum transport shows that rotating 
models including meridional circulation and shear turbulence 
are unable to reproduce the near uniformity of the solar ro- 
tation profile (Pinsonneaultetal. 19891 Chabover et al.1 119951: 



Eggenberger et_aL 20051; Turck-Chieze et alJ 20101) . The helio- 
seismic results stimulated various attempts to obtain similar ob- 
servations for stars with different masses and at various evolu- 
tionary stages. In the past few years, the spectrographs devel- 
oped for extra-solar planet searches have enabled the detection 
of s olar-like oscillations for a handful of solar-type stars (see 
e.g. Bedding & Kieldsen 2008), as well as for a few red giants 
( Fran dsen et al.ll2002l:lBarban et al.ll2004 120071: iDe Ridder etaD 
2006). Owing to the limited duration of these ground-based as- 
teroseismic campaigns, no clear additional constraints have how- 
ever been obtained about the internal transport of angular mo- 
mentum, with only one tentative evidence of rotational spl ittings 
being reported for the F9 V s tar ft Virginis dCarrier et al.| [2005: 
Eggenberger & Carrier 2006). We recall that rotation lifts the 
degeneracy in the azimuthal order m of non-radial oscillation 
modes. This results in (2d +1) frequency peaks in the power 
spectrum for each mode. The frequency spacings between these 
peaks are called rotational splittings and are directly related to 
the angular velocity and the properties of the modes in their 
propagation regions inside the star. 

Wit h the launch of the CoRoT dBaglin et alJ 120061) and 
Kepler dBorucki et alJl2010h space missions, solar-like oscilla- 
tions are now being detected and characterized for a very large 
number of stars. Photometric measurements of solar-like oscil- 
lations with the Kepler spacecraft have recently led to the pre- 
cise determination o f the rotational fre quency splittings of mixed 
modes in red giants dBeck et alJl2012h . Mixed modes are identi- 
fied in the power spectrum as dense clusters of modes referred 
to as 'dipole forests' for i = 1 modes. The mode at the cen- 
tre of such a forest is dominated by its acoustic character and is 
more sensitive to the external layers, while adjacent modes in the 
'wings' of the forest are more gravity-dominated and sensitive 
to the stellar core. In this work, we inv estigate whi c h con straints 
can be brought by the observations of Bec k et alJ (1201 2l) to the 
theoretical modelling of the physical processes transporting an- 
gular momentum in stellar interiors. 

The comparison between theoretical rotational splittings ob- 
tained with models including shellular rotation on ly and the 
splitti ngs observed for the red giant KIC 8366239 (Bec k et ail 
l2012h is first discussed in Sect. 2. Constraints imposed by these 
observed rotational splittings on an additional mechanism for the 
transport of angular momentum in stellar interiors are studied in 
Sect. 3, while our conclusions are given in Sect. 4. 



2. Models including shellular rotation only 

We use the Geneva stellar evol ution code that includes a detailed 
treatment of shellular rotation ( Eggenb erger et al.l2008l) . The ba- 
sic physical ingredients of numerical models of ro tating stars are 
described in Sect. 2 of Eggenberger e t alJ d2010bl) . We simply 
recall here that the transport of angular momentum obeys the 
equation 

where r is a characteristic radius, p the mean density on an iso- 
bar, O(r) the mean angular velocity at level r, and U (r) the ver- 
tical component of the meridional circulation velocity. The first 



term on the right-hand side of Eq. (fT) describes the advection 
of angular momentum by meridional circulation. In the classi- 
cal context of shellular rotation, both meridional circulation and 
shear mixing are considered as the main mixing mechanisms in 
radiative zones. The diffusion coefficient D then only accounts 
for the transport of angular momentum by shear instability with 
D = £> S hear- In convective zones, solid-body rotation is assumed 
to occur, in accordance with the solar case. 

A modelling of the red giant KIC 8366239 is performed 
by computing rotating models including shellular rotation only. 
These m odels are evolved with the solar chemical composition 
given by iGrevesse & Noelsl d 1993b and a solar calibrated value 
for the mixing-length parameter. We first adopt an initial veloc- 
ity of 50 km s _1 on the zero-age main sequence (ZAMS) and ob- 
tain a model with an initial mass of 1 .5 M in the H-shell burn- 
ing phase that reproduces the global asteroseismic properties of 
KIC 8366239 (i.e the large separation and the frequency of max- 
imum oscillation power). This model is characterised by a mass 
of 1.4993 M G , a luminosity of 12.59 L , a radius of 5.27 R and 
a mass of the inert helium core of 0.198M G . The rotation pro- 
file of this model is shown in Fig.Q] As expected from previous 
computations of rotating models including shellular rotation, the 
rotation profile during the red giant phase is very steep as a re- 
sult of the central contraction occurring at the end of th e main 
sequence dPalacios et al. 2006; Eggenberge r et al.ll20 10c). 

The rotational splittings of the dipole modes (i — 1) are then 
comp uted for this model using the Aarhus adiabatic pulsation 
code dChristensen-Dals g aard 2008). A mean value of 32.8//Hz 
is then found for modes at the centre of the dipole forests. This 
is more than three orders of magnitude higher than the mean 
value of 0.135 + 0.008 pHz observed in KIC 8366239. This huge 
discrepancy between observed and theoretical mean values of 
rotational splittings mainly comes from the large increase in the 
rotational velocity in the core of the model. Moreover, the ratio 
of the average rotational splitting of the first modes in the wings 
of the dipole forest to the mean splitting of the forest centre is 
equal to 2.5, which is significantly larger than the observed value 
of 1.5. We recall here that all mixed I = 1 modes are sensitive 
to the properties of the stellar core, but that they contain a differ- 
ent amount of pressure and gravity-mode influence. The mode 
at the centre of a dipole forest is most efficiently trapped in the 
acoustic cavity and is thus more sensitive to the properties of the 
stellar envelope than the other modes, which exhibit larger grav- 
ity components towards the wings of the dipole forest. The ratio 
of the rotational splittings of modes in the wings to those for 
modes at the centre of a dipole forest therefore gives important 
constraints on the difference between the rotation rates of the 
stellar core and the outer layers. The discrepancy between the- 
oretical and observed values of this ratio indicates that models 
including shellular rotation only, predict rotation profiles during 
the red giant phase that are too steep to correctly reproduce the 
measurements of rotational splittings in KIC 8366239. We must 
however emphasize that the computation of rotational splittings 
is done here by using a standard first-order perturbation expres- 
sion. While the assumption of slow rotation is valid in the enve- 
lope of our model, it is questionable in the very central layers, 
where the angular velocities become comparable to the frequen- 
cies of the oscillation modes. One can then wonder whether the 
large difference found between theoretical and observed rota- 
tional splittings is really related to the too-steep rotation profiles 
predicted by models including shellular rotation only, and is not 
simply due to an improper computation of these splittings for a 
model with a rapidly rotating core. 
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Fig. 1. Rotation profiles for models of KIC 8366239 computed 
with shellular rotation only. The continuous and dashed lines 
correspond to an initial velocity on the ZAMS of 50kms -1 and 
1 km s . The vertical dotted and dashed-dotted lines indicate the 
locations of the border of the helium core and the convective en- 
velope, respectively. 



To investigate this point, a model of KIC 8366239 is com- 
puted with a low initial velocity on the ZAMS of 1 kms" 1 . The 
corresponding rotation profile is shown in Fig.Q] A decrease in 
the initial velocity leads to a global decrease in the internal rota- 
tion rates. This results in a strong decrease in the rotational split- 
tings for modes at the centre of the dipole forests with a mean 
value of 1.6/iHz. This is still however more than one order of 
magnitude larger than the observed value. Moreover, the ratio of 
the splittings of modes in the wings to those for modes at the cen- 
tre of dipole forests remains significantly larger than observed 
(2.4 instead of 1.5). The angular rotational velocities in the cen- 
tral parts of this model are much lower than the modes frequen- 
cies so that the expression used for the computation of rotational 
splittings is perfeclty valid in this case. A slow initial rotation is 
thus found to produce rotation profiles during the red giant phase 
that are also too steep compared to observations. We conclude 
that meridional circulation and shear instability alone produce 
an insufficient coupling to account for the rotational splittings 
observed in the red giant KIC 8366239, and that an additional 
mechanism for the transport of angular momentun must operate 
in stellar interiors during the post-main sequence evolution. 



3. Constraints on an additional mechanism for 
internal angular momentum transport 

We now investigate in more detail how the rotational splittings 
observed in KIC 8366239 can constrain the needed additional 
mechanism for the transport of angular momentum. For this pur- 
pose, we introduce a new viscosity v a dd corresponding to an ad- 
ditional physical process for the transport of angular momentum 
in radiative interiors. The total diffusion coefficient in Eq. ([TJ is 
then given by D — D s h ear + ^add- The additional viscosity is not 
introduced into the equation describing the transport of chem- 
ical elements. This corresponds well to the two main mecha- 
nisms currently proposed to efficiently extract angular momen- 
tum from the central core of a solar-type star - magnetic fields 
and internal gravity waves - having a strong impact on the trans- 
port of angular momentum and only an indirect effect on the 
transport of chemical elements through the change of the rota- 
tion profile. As a first step to progress in our understanding of 
this additional physical process, we simply assume that v a dd is 
constant. We then investigate how its value can be constrained by 
asteroseismic observations of rotational splittings in the red gi- 



0.6 





2.5 



a 1.5 



1 T 
■W 2xl0 4 ; V,: 2Q 
".dd : 3x10*; V,: 20 



5xl0 4 ; V,: 20 

3xl0 4 ; V,: 40" 
5x10"; V,: 40 



H 1 1 1 1 1 1 h 




185 



190 

v [/uHz] 



195 



Fig. 2. Rotational splittings 5v for modes of the dipole forest 
near 190 //Hz for models of KIC 8366239 computed with differ- 
ent initial velocities on the ZAMS (Vj in kms -1 ) and viscosities 
(vadd in cm 2 s -1 ) corresponding to an additional mechanism for 
the transport of angular momentum in radiative zones. The bot- 
tom panel shows the rotational splittings normalised to the split- 
ting of the mode at the centre of the dipole forest. The dotted 
horizontal lines indicate the observed values of both the mean 
rotational splittings of modes at the centre of the dipole forests 
(top pannel) and the ratio of the splittings of the first modes in 
the wings of the dipole forest to the splitting of the mode at the 
centre of the forest (bottom panel). 
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Fig. 3. Rotation profiles in the central layers for the models 
shown in Fig. [2] 



ant KIC 8366239 and thus determine the efficiency of this mech- 
anism. 

A small grid of rotating models including this additional 
mechanism are then computed for different initial velocities on 
the ZAMS (Vi) and values of the viscosity (v a dd)- To correctly 
reproduce the rotational splittings observed in KIC 8366239, a 
viscosity v a dd = 3 x 10 4 cm 2 s _1 is found with an initial veloc- 
ity on of 20kms~ 1 . Figure [2] shows the comparison between the 
theoretical rotational splittings for this model and the observed 
values obtained for KIC 8366239, while the corresponding rota- 
tion profile is shown in Fig. [3] We note that only the splittings 
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near 190/vHz are shown in Fig.[2l but that the results are identi- 
cal for the splittings observed near 160, 175, and 200 //Hz. 

We find that the two input parameters of our models (V\ 
and v a dd) can be precisely determined thanks to the simultane- 
ous constraints coming from the observed values of the mean 
rotational splitting and from the ratio of splittings for modes in 
the wings to those for modes at the centre of the dipole forest. 
The bottom panel of Fig. [2] shows that the aforementioned ratio 
enables a direct and precise determination of v a dd- As shown in 
Fig. [3] an increase in v a dd results in a more efficient transport of 
angular momentum, hence a flatter rotation profile in the radia- 
tive zone. Since the ratio of splittings for modes in the wings to 
those at the centre of the dipole forest is sensitive to the dif- 
ference between the rotation rates of the stellar core and the 
outer layers, it is found to significantly decrease when the vis- 
cosity increases. The observed values of 1.5 for the first modes 
in the dipole forest is thus found to strongly constrain v a dd to 
3 x 10 4 cm 2 s _1 , independently of the initial velocity adopted in 
our models. The initial velocity V; of the model can then be ob- 
tained thanks to the value of the rotational splitting of the modes 
at the centre of the dipole forest. For a given viscosity, this value 
increases with the initial velocity as a result of the global in- 
crease in the rotation rates within the stellar interior. As shown 
in the top panel of Fig. [2] the observed rotational splitting of 
0.135 + 0.008 /iHz for modes at the centre of a dipole forest gives 
an initial velocity on the ZAMS of 20 km s _I (for the viscosity of 
3 x 10 4 cm 2 s , which is determined from the ratio of splittings 
for modes in the wings to those for modes at the centre of the 
dipole forest). The value of 20km s _1 found for the initial ve- 
locity on the ZAMS depends of course on the assumption made 
about the surface braking during the main sequence. A stronger 
braking results in a larger loss of angular momentum during the 
main sequence so that a higher value of the initial velocity is 
then required to correctly reproduce the observed mean value of 
rotational splittings. We note that the value obtained for the vis- 
cosity is insensitive to the adopted braking, because it is directly 
contrained by the ratio of splittings for modes in the wings to 
those for modes at the centre of the dipole forest, independently 
of the mean value of the splittings. We however recall here that 
the value of 3 x 10 4 cm 2 s _1 found for the viscosity in the radia- 
tive interior has been obtained by using the usual assumption of 
solid-body rotation in convective zones. 

We conclude that the efficiency of the additional mecha- 
nism for the transport of angular momentum in radiative zones 
operating during the post-main sequence evolution, can be 
strongly constrained by the observation of rotational splittings 
for modes in the wings of the dipole forests for the red giant 
KIC 8366239. The comparison between the determined viscos- 
ity of 3 x 10 4 cm 2 s and the coefficient D s h e ar shows that sim- 
ilar values are found for Z) s h ear at the border of the convective 
envelope, but that this coefficient rapidly decreases towards the 
centre of the star with typical values of lower than 50 cm 2 s _1 
in the stellar core. This underlines the need to introduce a new 
transport mechanism to efficiently transport angular momentum 
in the central parts of the star to obtain a core rotating at a moder- 
ate rate during the red giant phase. Interestingly, the viscosity of 
3 x 10 4 cm 2 s _1 deduced from asteroseismic observations of a red 
giant is similar to the value required to correctly account for the 
observed spin-down of slow rotat ing solar-type stars duri ng the 
beginning of the main sequence dDenissenkov et aLC OlO). This 
suggests that the same unknown physical process for the inter- 
nal transport of angular momentum may be operating during the 
main-sequence and the post-main sequence evolution. 



4. Conclusion 

We have found that models of red giant stars including shel- 
lular rotation only predict steep rotation profiles in the central 
stellar layers, which are incompatible with observational mea- 
surements of rotational splittings in the red giant KIC 8366239. 
Asteroseismic observations indeed show that meridional circula- 
tion and shear instability alone produce insufficient internal cou- 
pling and that an additional mechanism for the internal transport 
of angular momentum is thus needed during the post-main se- 
quence evolution. 

We then investigated in more detail how the rotational split- 
tings detected in KIC 8366239 can be used to characterize this 
physical process by computing stellar models including an ad- 
ditional viscosity related to this mechanism. We found that the 
ratio of splittings for modes in the wings to those for modes at 
the centre of the dipole forests strongly constrains the efficiency 
of this additional physical process and leads to a precise deter- 
mination of the corresponding viscosity of 3 x 10 4 cm 2 s . Such 
a value of viscosity may indicate that the same unknown physi- 
cal process for internal angular momentum transport is at work 
during the main-sequence and the red giant phase. 
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